A MODEL FOR QUARK-GLUON PLASMA WITH PENTAQUARK BARYONS 

AND TETRAQUARK MESONS 



R. S. Kaushal 

Department of Physics, Ramjas College (University Enclave), 
University of Delhi, Delhi-110007, India 
and 

Department of Physics & Astrophysics, University of Delhi, 
Delhi-110007, Indifl 

With a view to exploring a new kind of phase transition in the process of hadronization of quark- 
giuon plasma (QGP) we investigate the occurrence of pentaquark baryons and tetraquark mesons 
in the system. For this purpose, the frame work of an analoguous Saha's ionization formula for 
the colored ions in the system is used. The study of color-ionic-fraction (GIF) of multiply (color) 
ionized to unionized quark clusters (termed as "quarkons") as a function of temperature is carried 
. out. It is pointed out that not only the temperature of the fire-ball in the relativistic heavy ion 

collisions evolves with respect to space and time but also the GIF associated with a particular stage 
of ionization. Further, for the case of single color-ionization a correspondence of the present results 
with those available for the bubble nucleation mechanism in QGP is demonstrated. 
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I. INTRODUCTION 

It is now well known that heavy ion collision experiments performed at ultra-relativisitc energies 
(> 200 GeV/nucleon) can offer the same conditions as they existed at the time of origin of early 
universe, i.e., the conditions immediately after the big-bang. As a matter of fact the universe, during 
its space time evolution, has passed through a stage (though for a short period) in which a deconfined 
state of quarks and gluons (in brief called quark-gluon plasma or QGP characterized as a 'fire-ball') 
existed and the same is now expected to be formed during mini-bang, i.e., in the relativistic heavy 
ion collision (RHIC) experiments. In recent years, these studies have been the subject of great 
interest [1,2]. Particularly, the efforts have been there to look for the underlying mechanism and the 
processes which have led to the formation of hadronic state of matter out of the deconfined state or 
the vice versa. Various models to this effect have been discussed in the literature [3]. Until recently 
the picture considered was as follows: It is believed that initially existing deconfined state of quarks 
and gluons cooled down and as it evolved with space and time, the formation of hadrons (tri-quark 
baryons and two-quark mesons) took place via several stages of phase transitions. In particular, this 
quark-gluon state first converted into quark gluon soup[4,5], then into the diquark gluon plasma[6] 
(since diquark is also a colored object) and finally into the hadronic phase. 

It may be mentioned that beside the baryons (qqq) and mesons (qq) the existence of other (likely 
to be) stable quark clusters (now onward termed as 'quarkons') has been noticed only recently in 
various experiments [7]. Not only this, several attempts have already been made[8,9] to understand 
them in theoretical terms. While the theoretical understanding of these newly discovered objects 
(which in particular, are the pentaquark baryons like ^+(1540), C°(3099) and S (1860) and some 
of them are designated as (ududs), (ududc), (dsdsu) etc. and tetraquark mesons like (udud)) has 
yet to be established, their role in the physics of QGP has already become[10] of great interst. As 
a matter of fact the possibility of production of pentaquark 9^- baryon in the RHIC experiments 
has recently been discussed by Chen et al.[10] in a kinetic model via the processes KN*r^ 9, KN <-> 
7r6',and tiN ^K9. 

In fact once the existence of pentaquark baryons and tetraquark mesons is confirmed, the studies 
of QGP would require another stage of phase transition. As a result a new picture of QGP formation 
in RHIC experiments would emerge. The purpose of the present paper is to discuss one such model 
of QGP formation, of course, within the frame work of somewhat less refined scheme, i.e., within the 
frame work of color-quark chemistry which once played[ll] an important role in understanding the 
quark dynamics of hadrons. For this purpose, we shall explore here (perhaps for the first time to the 
best of our knowledge) the tools of Saha's ionization formula[12] (SIF) but now for color-charges in 
stead of the customary Coulomb charges. In particular, a computation of color ionic fraction (GIF) 
will be carried out for different quarkons and 'diquarkons' (diquark clusters). 

In the next section, we review the knowledge of SIF for coulombic ions and discuss the possible 
extension of results to the case of colored ions. In Section 3, we highlight some typical processes 
responsible for the occurrence of pentaquark and tetraquark systems and carry out the computation 
of Saha's color-ionic fraction. A possible connection between the single color-ionization process and 
the well-studied nucleation mechanism for the hadronic phase transitions, is highlighted, perhaps for 
the first time, in Section 4. Finally concluding remarks are made in Section 5. 

II. BRIEF REVIEW OF SAHA (COULOMB) IONIZATION FORMULA (SIF) 

A detailed study of ionization of gases by thermal excitation was carried out by M. N. Saha about 
eighty years ago[12]. Particularly, the question as to what happens when a gaseous mass consisting 



3 



of atoms is heated to a very high temperature was discussed by Saha and a formula for computing 
the fraction of the number of (r + l)-times ionized to the r-times ionized atoms was derived in terms 
of pressure, temperature and internal energy of the system. 

The generalized version of the SIF for the chemical reaction J2r<^r^r ^ Ss^s-^s, as an outcome 
of classical statistics and the law of mass action provides, (cf. Ref.(12), p. 658) 



lnKp^—— + -lnT + lnI + lnZ, (1) 
Jril 2 

where the fractions, -ftTp, /, and Z of corresponding partial pressures (P's), intrinsic constants (I's), 
and parition functions (Z's) for products and reactants in the given chemical reaction, are defined 
as 



Here U and T respectively are the internal energy and the temperature of the system and R is 
the gas constant. Further, one restricts to the electronic excitations of the atoms. In the present 
context however we shall consider, in analogy, the quark excitations of the quarkons and deal with 
the transfer of one unit of color (anticolor) charge at a time from parent to daughter quarkon or the 
vice- versa. 

It may be noted that for a simple atomic process for an element M, the ionization and capture 
described by M, M # M+-|-e~, (M+ is the ion and e~ is the electron), equation (1) is written as[12] 



where is the electron spin multiplicity,i.e.,2 and Z^ denotes the partition function of the atom in 
the argument , and the constant n is given by k = /n((27rm)^/^A;^/^//i^). Under the assumption that 
one initially starts with the element M and then goes on heating the system in a confined space 
to a temperature T, the fraction x of the ionized atoms, also linked with the partial pressures Pm, 
Pm+jPc- through Pg = Pm+ — nxkT, Pm — n{l — x)kT, can be computed from (3) as 



where P = P^- +Pm+ +Pm is the total pressure and the contribution of the last term in (3) is ignored. 
Further, formula (3) can also be written for the ionization of the atom at any stage of ionization, 
namely from to M^~^^ (where denotes the atom which has already lost r-electrons) , as 



In Pe = —-— + -InT + K + ln[2— ————). (5) 

Here Ur is the heat of ionization from to M'^'^^ and the electronic partition function of M^", 
Ze(M''), for the r-times ionized atom is given by 



s 



(6) 
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where gr and grs are being the weights corresponding to the ground and s— th excited state of M**; 
Xs is the excitation energy of state s, and the temperature T is now onward expressed in the units 
oik. 

For the sake of ready use we give below another version of formula (1), namely[13] for the reaction 
— 0, one writes (1) in the form 

K{P, T) = P- ''Ke- '''^'1'^ (7) 

which, for a paricular reaction + e~" — A° = with va+ = 1, I'e- = 1; = — 1, reduces to the 
form 



or after using P — {N/V)T, one writes (8) as 

K{P^T) = ^.M±(^)3/2_y3/2g-X./T_ (9) 

In this case, if one defines the degree of ionization (ionic fraction) a — N^/Nq, then eq.(9) can be 
recast (see, Rumer and Ryvkin, Ref.(13)) as (jz^ = 20(7", y), leading to 



a(T, V) = -0(r, V) + (0(r, V)' + 20(r, v))'/' (lo) 

where 

cl>{T,V) = {T/Tor/\e-^^/^,withTo = (-^)(2^^)^/^ (11) 

27rm geg+ V 

is a positive, dimensionless function. Note the difference between the definitions of fractions a in 
eq.(lO) and x of eq.(4). We shall however investigate a as a function of T. 

We shall restrict here only to two-body final states and use ge — ^ for quarks and ge = ^ for scalar 
diquarks. Note that the formula (10) is more convenient if the ionized gas is confined to a fixed 
volume as is the case with fire-ball in RHIC experiments. Further, we shall use the version of SIF as 
described in eqs. (8)-(ll) for quarkons and diquarkons, which in certain situations will correspond 
to pentaquark baryons and tetraquark mesons. The picture considered here is as follows: 

Whether it is deconfined phase of quarks and gluons that existed at the time of early Universe or 
the same is produced in RHIC experiments in a localized region (fire-ball), the system cools down 
and evolves in space and time thereafter. As a result, the hadronization takes place after the system 
passes through (may be for a short while) a diquark-gluon phase or a mixed quark-diquark-gluon 
phase. Alternatively, one can think of (at least in the case of RHIC experiments in the forward 
process) the breaking of the normal hadrons at such high energies into their constituent quarks and 
thereby simultaneously forming a metastable state of quark and/or diquark clusters. We assume that 
quarks or diquarks in these clusters are bound by Coulomb-hke color forces, as the corresponding 
potential is found to work well for quark constituents of the nucleon[14] and also of pentaquark 
baryons[8]. Note that pentaquark baryons or tetraquark mesons could be one of the possibilities of 
these quarkons or diquarkons present initially. In the next stage , these clusters will decay and give 
rise to a mixture of quarks and diquarks and finally to quarks only. With regard to this picture. 
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we make the following simplifying but plausible assumptions: (i) Role of gluons in this model will 
appear only through the transfer of color charge; (ii) The transfer of one quark will carry only one 
unit of color charge ( as is the case with the Coulombic charge of electron) and that of a diquark 
will carry two units of color charge; (iii) No account will be made of fractional Coulomb-charge on a 
quark or diquark (or on the quarkon or diquarkon for that matter); (iv) We also ignore an account 
of nature of color on a particular quark and also that of its flavor for the time being, (v) While 
neglecting other production or recombination processes, we shall concentrate here only on the decay 
of quarkons and diquarkons into lighter quarkon channels -perhaps more justified for the formation 
of QGP. Some of these assumptions conform to the spirit of color-quark-chemistry studied earlier by 
Chan Hong- Mo et al.[ll]. 



III. QUARKONS AND DIQUARKONS IN THE FIRE-BALL AND THEIR DECAY 
A. The processes involving pentaquark baryons and tetraquark mesons 

The quarkons and diquarkons, once formed in the fire-ball in RHIC experiments, subsequently 
decay into quarks and diquarks and in due course contribute to the deconfined state of these objects 
and that too in a confined volume. While the studies can be easily extended to a general case , 
for simplicity we concentrate here only on the decay of pentaquarkons and tri-diquarkons, denoted 
respectively by (5(o) ^^'^ -^(o)' Here superscript denotes the number of quarks(diquarks) present 
initially in the cluster (termed as 'parent') and subscript represents the number of quarks (diquarks) 
left over after the decay (termed as 'daughter'). In other words, the quarks (or for that matter 
diquarks) released finally from the cluster due to the increase in temperature of the fire-ball, will be 
able to keep a track of their parentage. The possible processes associated with the decay of these 
systems can be listed as follows: 



V(0) ^ V(4) + 



Q(3) + Q{2) ~^ Q(3) + 2(5(1), 



(12) 
(13) 



Q[2) + Q(3) ~^ Q(2) + 3(5(1), 



'(3) 



(2) 



'(1)' 



2(5(2) Q{1) ~^ ^Q(l): 



(1)' 



V(0) ^ V(3) ^ ^{1)1 



(14) 
(15) 
(16) 



(5(2) + 2(5(1), 



2(5(2) 



(17) 
(18) 



(19) 
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1 


^(0) 






^(0) 


^(2) + 


^(3) 
^(1) 



(20) 
(21) 

(22) 

-3L>g, (23) 

dU - 22?a)> (24) 

-D[o) = Q[o) ~^ (25) 

While the processes (22)- (24) will contribute to the diquark or quark-di quark plasma, the processes 
(12) to (21) and (25) will be responsible exclusively for quark plasma. In fact the role of pentaquark 
baryons or tetraquark mesons, if at all manifests in the formation of QGP, it is reasonable to assume 
that the same will appear through the decays (color ionization) of pentaquarkons <5(o)' tetraquarkons 

Q(o)' '5(4)5 -D[o) including that of other parent clusters; whereas the decay of triquarkons Q^^-^, Q(3), 

diquarkons Q(2)' Q(2)' monoquarkons Q^^^, Q(i) will contribute to quark and diquark 

plasma. We compute here the color ionic fraction (GIF) for some representative cases using SIF for 
the colored ions and study the same as a function of temperature at various quark number densities. 

B. Calculation of color ionic fraction for various processes 

When computing the GIF using (10) and (11), we need some ingredients about the quark compos- 
ites. The same we use from our earlier works[8,14] in which a quark diquark (QDQ) model for the 
nucleon[14] and a quark double diquark (QDDQ) model for the pentaquark baryons[8] are proposed. 
With regard to the geometry of the fire-ball (though it depends on the colliding ions) in a typical 
RHIG experiment, the results are taken from Karsch and Petronzio [15] and others[16] derived in 
the context of J/ip suppression. 

In particular, we concentrate on some typical processes (12), (13), (18), (19) and (21). In fact 
these arc the cases for which some parameters about the concerned quarkons/diquarkons are 
already known[8,14]. Note that in the same spirit we consider here Goulomb-like color forces 
among the quarks in a quarkon and among the diquarks in a diquarkon and use the hydrogenic 
model for deriving the various parameters from the experimental results. In this way the ion- 
ization potential,Xs, in (11) for various clusters are computed. Thus, we use the following ingredients: 

QDQModel[14] : 

rUg = blS.OMeV; mo = GSl.SMeV; 6oc = 0.47/m; 6o^, = 0.57/m; 

le^l = 349.0Me1/; |e^| = 215.8Me1/; pc = 1-64; py = 1.23. (26) 

QDDQModeI[8] : 
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= 0.8/m; So„ = 0.48/m; 5^ = 0.723; 5^ = 1.13; 

= 90.4Mey; \El\^ 235 AMeV. (27) 

For the tetraquark (a state of diquark-antidiquark (DD) bound system) case, we assume the 
following relationships for couplings: 

f^qq _ ^DD 
fiqq SdD ' 

and derive Sj^j^ using the models described in Refs.(8) and (14). This leads to Sj^i) = 1.06 for 
light quarks and Sj:,i)_ = 0.56 for the heavy (charm) quarks and accordingly for the tetraquark 
meson we find[8] \E^^\ = IQlAMeV, for the hght DD-case, and \E^^\ = 54.2Mey,for the heavy 
(charm)£)£)-case. Next we calculate the GIF, a, from (10) for the following four cases: 

Casela: Decay of the pentaquarkon (cf. eq.(12), valence) :(5(oj ~^ ^(4) + Q(i) 

In this case the quark is released from the valence in the QDDQ model of the pentaquark 
baryons[8] and the process will contribute to the quark gas (cf. Fig. la). 

Caselb: Decay of the pentaquarkon (cf. eq.(13), core):(5|o) Qfs) + Q(2) 

This is the case in which one of the diquarks from the core in the QDDQ model[8] is released. 
This process will contribute to the diquark gas (cf. Fig. lb). 

Casella: Decay of the triquarkon (cf. eq.(19), valence) :Q|q| Q^^^ + Q^^^ 

In this case the triquarkon is ionized and will contribute to both quark and diquark gases or to 
their mixture. The valence in the QDQ model of nucleon[14] is released (cf. Fig. 2a). 

Casellb: Decay of the triquarkon (cf. eq.(19), C0Te):Q^^^ — > Q(2) + 

This is the case in which the quark from the diquark core of the nuclcon[14] is released and will 
contribute to both quark and diquark gases or to their mixture (cf. Fig. 2b). 

Caselll: Decay of tetraquarkon (cf. eq.(18):Q^^^ — > Q^) + Q(2) 

In this case, tetraquarkon breaks up into two diquarks and will exclusively contribute to diquark 
gas and the model of Ref.(14) is used (cf. Fig. 3). 

CaseIV:Decay of diquarkon (cf. eq.(21):(5(o) ^ Q(i) + Q(i) 

This process will exclusively contribute to quark gas and the parameters are derived by assuming 
the diquark to be the same as in the core of the nucleon in the QDQ model[14]. 

Corresponding to these processes the GIF is computed from eq. (10). The values of various 
parameters used in the calculations are those given in eqs. (26) and (27) and in Table 1 for a sample 
case of O^^-O^^ coUision. For this case, an order of magnitude estimate for the size of the fire-ball 
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TABLE I: Values of various quantities used in the calculations of color-ionic fraction, a, from eqns.(lO) and (11). 



process 




9+ 


9 a 


m(Afel/) 


XsiMeV) 


To{MeV) 


Case la 


2 


1 


2 


400.0 


235.4 


3026.0 


Case lb 


1 


2 


2 


460.0 


90.4 


1664.0 


Case Ila 


2 


1 


2 


288.3 


215.8 


4225.0 


Case lib 


2 


1 


2 


333.3 


349.0 


3656.0 


Case III 


1 


1 


1 


340.5 


191.4 


2257.0 


Case IV 


2 


2 


1 


250.5 


349.0 


1935.0 



and of (Nq/V) are carried out using the parameters given in Refs.(15) and (16). The calculated 
reduced mass, m, of the two-body final state and the spin multiplicities ge, g+, Qa along with the 
calculated ionization energy Xs,are listed in Table 1 for different cases. Further, for the above six 
cases the calculated results for the GIF, a, as a function of temperature T are shown in Figs. 1-4 
corresponding to the three typical values of the particle number density {Nq/V)^ namely 10x10^ 
(continuous curve), 50x10^ (dashed curve) and 90x10^ MeV^ (dotted curve). While T = 2QQMeV is 
believed to be the temperature for the hadronization, the a is shown in these figures in the range 
T = 50MeK to 350Mey. Some crucial difference around this temperature can be seen from these 
figures for differnt cases and for different particle number densities. 



IV. A CONNECTION BETWEEN COLOR-IONIZATION AND NUCLEATION MECHANISM 

Nuclcation or precipitation mechanism in a medium has been the subject of study in various field 
and in different contexts for many decades now[17-19]. The nucleation rate in the context of QCD and 
QGP has been investigated by Csernai and Kapusta[20] among others[21]. In this case one computes 
the probability that a bubble or a droplet of 74-phase appears in a system initially in the S-phase 
near the critical temperature. Again for the case like early universe or RHIC studies homogeneous 
nucleation theory is found more convenient. In fact a droplet of critical size is metastable, it is 
balanced between evaporation and accretion. Somewhat similar is the picture in the present case 
in which a quarkon (or a diquarkon) of critical size ( in the sense that it becomes a normal hadron 
for a certain number of quarks(diquarks) and a part of the medium just for another number) plays 
the role of a droplet or bubble. The decay of quarkon or diquarkon in this case, can be considered 
as the process of evaporation of a droplet. Further, for the study of nucleation rate both classical 
theory[17,18], based on thermodynamics, and the modern theory of Langer[19], based on statistical 
methods have been used. The concept of nucleation from one vacuum to another has also been 
extended to the domain of quantum field theory by several authors[22]. 

The interesting part in all these studies in rather disconnected fields is that the derived expression 
for the nucleation rate, in general, has a common feature in terms of physical as well as mathematical 
contents. While the probabihty per unit time per unit volume to nucleate the dense liquid' phase 
from a dilute 'gas' phase is broadly expressible as 

I = /o.e-^=/^, (28) 

the prefactor /q (having dimensions of T^) here in general has a break-up into statistical and dy- 
namical parts in all the theories of nucleation. For example, in the Langer's theory[19], one writes 
the nucleation rate as 
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I = f f^o.e-^^/^, (29) 

where AF is the change in the free energy of the system due to the formation of critical droplet, 
and flo and k, are the statistical and dynamical prefactors which respectively are the measures of the 
available phase space volume and the exponential growth of the critical droplet. For the theories 
of early Universe while Jo in the literature( see, for example, [21]) is given by Iq — {Wc/2'kTY/'^T^ , 
Csernai and Kapusta within the frame work of a course-grained effective field theory approximation, 
obtain an expression for Iq in the QGP context as[20,21] 



where r]q = 14. 4T^, is the shear viscosity in the plasma phase;^^ is the correlation length in this 
phase; a is the surface free energy; is the radius of the critical-sized bubble, and Aw is the 
difference in the enthalpy densities of the two phases. 

It is worthwhile to compare the formula (29) (or (28)) with Iq given by (30) with the expression 
(9) defining the fraction of partial pressures for the products and the reactants, particularly for 
the single color-ionization. First of all note that beside the factor iV/N) on the right hand side 
of (9), the factor {gf..g^/gA has a statistical origin and corresponds to fig in (29). The remaining 
factor (27rmT//i^)^/^ is analoguous to k in (29) as it has roots in the dynamics. Secondly, a little 
simplification of the prefactor /q in (30) yields 



7o = (T/To)V^ withf, = i^^^^m^f'- (31) 
On the other hand, the prefactor, Kq, in (9) can also be expressed as 

Ko = (^/fo)'/^ (32) 

where Tq is given by 

r„ = (^f/=(^), (33) 

Vgeg+ 2-Km 

and K can still be shown to have dimensions of like / in (28), at least for the case of single 
color-ionization where K turns out to be a measure of pressure. Further it may be of interest to 
compare the magnitudes of the ionization energy Xs (cf. Table 1) and that of AF in (29). In fact 
they turn out to be of the same order (cf. Ref. (21)) in certain models. 

Thus, seemingly different mechanisms of single color-ionization and bubble nucleation have the 
same mathemathical content at their computational level. 



V. DISCUSSION AND CONCLUSIONS 



From the point of view of investigating the role of newly discovered pentaquark baryons and 
tetraquark mesons in the formation of QGP, we have made a modest attempt, perhaps for the first 
time, to demonstrate the viability of SIF for the case of colored ion systems. While the method 
is quite general for the study of multiply ionized systems, the case of single (color-) ionization is 
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investigated in detail. In spite of several assumtions made to simplify the computation of GIF, 
this latter quantity as a function of temperature yields the behaviour as expected. Interestingly, 
for this single-ionization case a connection of the prsent approach with the well studied problem 
of bubble nucleation in the context of QGP is demonstrated. There appears to be a one-to-one 
correspondence between the working formulae in the two approaches, in spite of different physical 
inputs in the two cases. To be more specific, our findings can be summarized as follows: 

(i) The fact that GIF, a, in all cases (cf. Sect. 3. 2) at different particle number densities ap- 
proaches to unity in the large-T limit, implies a complete dissociation of quarks or diquarks from 
the corresponding clusters, i.e., the formation of a noninteracting quark or diquark gas. 

(ii) Note from Figs, la and lb that at a given temperature, the contribution to diquark gas from 
the diquarks in the core of a pentaquarkon is more pronounced compare to the one obtained from 
its valence quark within the frame work of the QDDQ model. This conforms to the fact that the 
valence of antiquark in the pentaquark baryon is more tightly bound compare to the double diquark 
core -an important outcome of the QDDQ model[8]. 

(iii) In case of triquarkons, the core and valence (cf. Ref.(14)) contributions to both diquark and 
quark gases are comparable as far as the variation of a with T is concerned. However, the difference 
manifests in the low-T limit (cf. Figs. 2a and 2b). 

(iv) As expected, the diquarks of tetraquarkons contribute somewhat more to the diquark gas(cf. 
Fig. 3) compared to the diquarks of triquarkons (cf. point (iii) above). 

As the colored-ion system (such as the fire-ball in RHIC experiments) evolves in space and time, 
the temperature (cf. Refs.(15) and (16)) and hence the GIF (cf. eq.(lO)) also varies inside the fire- 
ball with space and time. This kind of space and/or time dependence of a can easily be investigated 
in the present frame work of SIF and that too for different stages of multiple ionization of quarkons 
and diquakons. However, for siplicity we have taken T as uniform all through out the volume of the 
fire-ball. 

Note that it is only the single (color-)ionization in the generalized version (1) of SIF which cor- 
responds to the bubble nucleation mechanism in QGP (cf. Refs. (20) and (21)). This shows the 
richness of SIF approach over the nucleation one as far as the hadronic-bubblc formation in QGP is 
concerned . Moreover, there is enough scope to investigate further higher-order phase transitions in 
QGP (if at all they exist) in the present approach by way of studying the two- or higher-stages of 
ionization of quarkons and/or diquarkons. Such studies are in progress. 
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FIG. la: Contribution to quaxk gas from pentaquarkons (valence): The results are shown for three values of paricle number 
densities (No/V), namely 10x10®, 50x10® and 90x10 Me\^^, respectively by continuous, dashed and dotted curves. 
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FIG. lb: Contribution to diquark gas from pentaquarkons (core). The description of curves is the same as in Fig. la. 
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FIG. 2a: Contribution to quark and diquark gas from triquarkons (valence). The description of curves is the same as in 

la. 
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FIG. 2b: Contribution to quark and diquark gas from triquarkons (core). The description of curves is the same as in Fig. 
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FIG. 3: Contribution to diquark gas from tetraquarkons. The description of curves is the same as in Fig. la. 
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FIG. 4: Contribution to quark gas from single diquarkons (core- like). The description of curves is the same as in Fig. la. 
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